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Abstract. MAGIC and H.E.S.S experiments are the first to produce images of 
supernova remnats (SNRs) in TeV gamma-rays. The gamma-radiation are produced 
either by electrons (due to inverse-Compton scatterings) or protons (due to pion 
decays) . We present a method to synthesize gamma-ray images of Sedov SNRs due 
to hadronic emission. The model is developed in the frame of a classic approach 
to proton acceleration and hydrodynamics of the shocks in a uniform interstellar 
medium; it includes energy losses of relativistic protons due to pp interactions. 
Our calculations show that these losses are important only for large densities of 
protons as it could be in case of interactions of the supernova shock with molecular 
cloud. Numerical simulations are used to synthesize radial profiles of hadronic TeV 
gamma-rays. 

Introduction 

Cosmic rays (CRs) are the widely studied. Supernova remnants (SNRs), main source of galactic 
CRs, are excellent objects to study magneto-hydrodynamics of nonrelativistic shocks and acceleration of 
cosmic rays, namely protons and electrons. These particles radiate from radio to 7-rays due to different 
types of emission. Experiments in high-energy astronomy observe all types of these emission. 

Most of galactic cosmic rays are belived to be produced by the forward shocks in SNRs. In particular, 
efficient proton acceleration changes the structure of the shock front and makes plasma more compressible 
that leads to lower adiabatic index, to increased shock compression factor and to some observed effects: 
reduced physical separation between the forward shock and the "contact discontinuity" (or reverse shock) 
[e.g Warren et at, 2006]; concave shape of the energy spectrum [e.g Reynolds and Ellison, 1992]; growth 
of some turbulence modes and to MF amplification in the pre-shock region [e.g Bell, 2004]; "blinking" 
X-ray spots originated from such growth of MF [Uchiyama et al, 2007]. 

Observations are expected to confirm that protons are accelerated in SNRs to very high energy and 
emit (TeV) 7-rays. Nevertheless, analysis of the broad-band spectra of SNRs shows that both electrons 
and protons may be responsible for TeV gamma-rays [e.g. SN \QQQ:Acero et al., 2010]. 

The properties of the thickness of the radial profiles of hard X-ray brightness are used to estimate the 
strength of the post-shock magnetic field [Berezhko et al., 2003]. Radial profiles of the radio brightness 
may constrain the time evolution of the electron injection efficiency [Petruk et al, 2011a]. In a simple 
fashion, the radial profiles of hadronic TeV 7-ray brightness are sensitive to the density of ambient 
medium. This property is the subject of the present study. 

Properties of the nonthermal images of Sedov SNRs due to radiation of accelerated electrons in radio. 
X-rays and 7-rays are systematically studied in [Reynolds, 1998, 2004] and [Petruk et al., 2009, 2011b, 
Papers I and II respectively]. Numerical models for synthesis of maps of adiabatic SNRs in uniform 
ISM and and uniform interstellar magnetic field (ISMF) from basic theoretical principles as well as their 
approximate analytical descriptions are developed in these papers. The main factors determining the 
azimuthal and radial variation of surface brightness of Sedov SNRs are determined there. 

These papers, as the present one, are limited to the test-particle approach because the non-linear 
theory of diffusive acceleration is not developed for shocks of different obliquity, while the obliquity 
dependence of various parameters is important for image modelling. 

In the present paper, we study properties of the radial profiles of surface brightness in 7-ray due to 
proton-proton interactions including of the energy losses of proton. 

Model 

Our model closely restores that used in Papers I and II. Let us consider an adiabatic SNR in 
uniform interstellar medium (ISM) and uniform interstellar magnetic field (ISMF) . We use quite accurate 
approximate formula in Lagrange coordinates [Petruk, 2000] for description of hydrodynamics of SNRs 
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in the adiabatic stage of evolution (Sedov solutions [Sedov, 1959]). Magnetic field is described following 
[Reynolds, 1998]. We do not consider amplification of the ambient field. 
At the shock, the spectrum of acceleration protons is taken as 

N,{Ep) = KsE-^ exp f-^r^) (1) 

where Ep^max is the maximum energy of protons, s is constant and we use s = 2. 

We assume that variation of the maximum energy with obliquity angle (angle between the ambient 
magnetic field and the shock velocity) is constant and the injection efficiency of protons is isotropic. 

The surface brightness is calculated integrating emissivities along the line of sight within SNR. 
Hadronic 7-rays appear as a consequence of the neutral pion and r^-mcson decays produced in inelastic 
collisions of accelerated protons with thermal protons downstream of the shock; the spatial distribution 
of the target protons is simply proportional to the plasma density. The hadronic 7-ray emissivity is 
calculated as [Kelner et al, 2006] 

1 

/dx 
app{E^/x)N{E^/x)Fj{x,Ej/x) — , (2) 



where x = E^/Ep, the cross-section is [Aharonian, Athoyan, 2000]: 

app{Ep) = 28.5 + 1.8 In {Ep/l GeV) mb, (3) 
the function is [Kelner et al, 2006] 

In (.'/;) / 1 — x^"' 



F^{x,Ep) = B^ 
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where 



= 1.30 + 0.14L + 0.011^2, (5) 

(6) 



1.79 + 0.11i + 0.008X2' 

1 

0.801 + 0.049L + 0.014L2 ■ 



and L = ln(i;p/l TeV). 

Energy losses of protons due to proton-proton interactions 

Modelling the surface brightness distribution and maps due to proton collisions is important to 
accurate the energy losses due to pion productions. The losses due to proton collisions are important for 
higher densities of target protons. It may be shown the proton collision losses may be described as 



dEp 
IT 



= 3KcnH(^ppiEp)Ep^kin, (8) 

pp 

where the factor 3 accounts for the production of 7r°, 7r+ and tt^ mesons, respectively and c, are 
speed of light, the proton target density, app^ET^, Ep) the differential cross-section for the interaction of 
two protons, Ep the energy of primary proton, 3k; = 0.51 if «; = 0.17 [Aharonian et al., 2000]. 

Let us compare the energy losses due to proton-proton interactions with the radiative losses of 
electrons. The losses of electrons are given by 

i..crjL)7i?). w 



'^^ J rad ^ \meC^ J \StT J ' 

where or is the Thomson cross-section, nie the mass of electron. The ratio of electron to proton radiative 
losses is -22 

Ep pp nnEp^TeV 
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where we used app w 33 mb, B^a is the magnetic field in lO"^^;, and Ee^xeV and Ep^TeV are the energy 
of electrons and protons in IQ^^eV . The losses of electrons end protons is similar when the dcinsity of 
protons 440 cm~^ and the typical galactic magnetic field Bi^g = 3 the maximum energy of electrons end 
protons are 30 and 100 TeV respectively. One can see that the losses of protons with energy 1000 TeV 
are comparable to losses of electrons with energy 30 TeV in magnetic field 3 fiG, if the number density 
of target protons is rather high 440 cm~^. 

Downstream evolution of the proton distribution 

Let the energy of proton at the time U, when it leave the region of acceleration, is Epi. Then it is 
smaller at the present time t, 

Ep = Epi£ad{aY^''^£pp{Ep, a), (11) 

because the terms responsible for the adiabatic E^^ and coUisional £pp losses are equal or smaller than 
unity; a = a/ a the Lagrangian coordinate, R the radius of SNR, 

£ad{o) = n{aY'\ (12) 

where n = n/us, index "s" denotes the value immediately post-shock, 

£pp{Ep, a) = {Ep/lGeVf-''^''^ I{a), (13) 

/i(o) and 1(a) are dimensionless self-similar functions 



li{a) = exp 



dx 



I{a) = exp 



c 



dx 



(14) 



(15) 



where C, = 5tei/2 = 1.21 • lO^^tan/zs, = f/lOOOyrs. It is clear from here that £pp is effective only 
where the density of target protons is large, at least urs ~ lO^cm"^. 

The energy spectrum of protons downstream of the shock evolves self-similarly 

-1 



Np{Ep,d,t) = K{d,t)E-' fi{d)£pp{Ep,dy 



X exp 



Ep,max£ad{aY^^^£pp{Ep, a) 



(16) 



with K{d,t) = KsK{d), K{d) = a3^/2fi(a)i+''(«)(«-i)/3. The downstream distribution of relativistic 
protons are modified by the adiabatic expansion of SNR. Losses due to inelastic collisions affects the 
distribution only when C, is not small, i.e. when uhs is large. 



Results 

Analysis of the formula (8) shows that the thickness of radial profiles of the TeV 7-ray surface 
brightness distributions are only function of density of target protons. If density of target protons is 
small then we can neglect the energy losses due to proton-proton interactions. However, when the 
density of target protons is more than ^ 400 cw? one can not neglect these losses. However,they reveal 
themselves in the radial profiles for even higher densities. 

In the fig.l, we show the influence of target protons density on the thickness of radial profiles of 
surface brightness. If the density increases the energy losses of protons increase as well and the thickness 
of the radial profile decreases. Since the downstream density is proportional to the pre-shock density, 
thus may be used as a probe of the ambient medium density. 



Conclusion 

We consider an influence of target protons density on the radial proflles of TeV 7-ray surface bright- 
ness. Downstream evolution of the proton distribution taking into account the energy losses due to 
proton-proton interactions is described. We show the influence of the target protons density on prop- 
erties of radial proflles of 7-ray brightness distribution. These properties can be used to estimate the 
density of molecular clouds that interact with the shock wave of SNR. 
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Figure 1. Radial profiles of the 7 -ray surface brightness due to hadronic emission for different densities 
of the target protons: lcm~^ (line 1), lO^'cm"^ (line 2), lO^cm"^ (line 3), 10'' cm"^ (line 4). Emax = 1000 
TeV, energy of gamma-rays is 1 TeV. 
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